Logarithmic Conformal Field Theory and Boundary Effects in the Dimer Model 
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We study the finite-size corrections of the dimer model onooxJV square lattice with two different 
boundary conditions: free and periodic. We find that the finite-size corrections in a crucial way 
depend on the parity of N; we also show that such unusual finite-size behavior can be fully explained 
in the framework of the c = — 2 logarithmic conformal field theory. 
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Universality, scaling, and exact finite-size corrections 
in critical systems have attracted much attention in re- 
cent decades [J, 0, H LI- It nas been found that crit- 
ical systems can be classified into different universality 
classes so that the systems in the same class have the 
same set of critical exponents, universal finite-size scaling 
functions, and amplitude ratios 0,0- Two-dimensional 
critical systems are parameterized by the central charge 
c Q, directly related to the finite-size corrections to the 
critical free energy 0. 

In this Letter, we address this question for dimers 
defined on a square lattice, with three main purposes: 
(i) we dissipate the confusion existing in the literature 
about the value of the central charge [7], |8| due to the 
(mis) interpretation of the finite-size corrections in terms 
of the central charge rather than the effective central 
charge; (ii) we give a bijection of dimer coverings with 
spanning tree and Abelian sandpile model, which not 
only allows a proper understanding of the dimer model 
but proves also very useful to calculate finer effects, like 
the change of boundary conditions; (iii) using this bijec- 
tion, we clarify and explain why a change of parity of the 
lattice size causes a change of the effective central charge 
but not of the central charge itself. 

The dimer problem on planar lattices belongs to the 
class of "free-fermion" models 0- Its solution has been 
obtained with the Pfaffian approac h [Toll and then re- 
produced by a variety of methods [lj. In contrast to 
the statistics of simple particles, the critical behavior of 
the dimer model is strongly influenced by the structure 
of the lattice space. The square lattice dimer model is 
critical with algebraic decay of correlators ^1 ■ F° r the 
dimer model on the anisotropic honeycomb lattice, which 
is equivalent to five-vertex model on the square lattice 
[l3| , the free energy exhibit a KDP-type singularity. For 
the triangular and some decorated lattices, the dimer 
model exhibit Ising-type transitions Q. Thus, it ap- 
pears that the dimer model itself has not a single critical 
behavior, but several critical behaviors associated with 
different universality classes. 

In what follows, we consider the finite-size effects for 



close-packed dimers on finite square lattices, with free 
boundary conditions on all sides (strip geometry), and 
with periodic boundary condition in one direction (cylin- 
der). In all cases we find them to be consistent with a 
central charge c = — 2. This conclusion relies on a careful 
distinction between the central charge c and the so-called 
effective central charge c e g — c — 24/i m ; n 01 > which is a 
boundary dependent quantity. The value c = — 2 is fur- 
ther confirmed by calculating the effect of a change of 
boundary conditions. Wc find that, in the scaling limit, 
it corresponds to the insertion of a boundary primary 
field of weight —1/8, belonging to a logarithmic confor- 
mal field theory with c = — 2. 

Finite-size analysis. - Let us consider the dimer model 
on an M x N square lattice C with M rows and N 
columns. The topology of C is fixed by the boundary 
conditions: it forms a rectangle if free boundary condi- 
tions are imposed in two directions, a cylinder or a torus 
if periodic boundary conditions are chosen in only hori- 
zontal direction, or two directions. 

The partition function of the dimer model is given by 



Zm,n(z v , Zh) — ^ 
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where the summation is over all dimer covering configu- 
rations, z v and Zh are the dimer weights in the vertical 
and horizontal directions, and n v and rih are the numbers 
of vertical and horizontal dimers. 

The partition functions of the dimer model with the 
boundary conditions discussed above can all be expressed 
in terms of Z a ,p(z, M, N) for a,/3 = 0A 4] with 
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Here z = Zh/z v , which we set equal to 1 from now on. 

The general theory about the asymptotic expansion of 
Z a p for large M, N has been presented in Q- For what 
follows, the asymptotic expansion of the free energy per 
unit length associated to Z a< p is all we need. The result 
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where z = l,z 2 = —2/3, ... ; G — 0.915965 is the Cata- 
lan constant, and the B p (a) are the Bernoulli polynomi- 
als, B 2 (a) — a 2 — a+\- These formulae allow to compute 
the asymptotic expansion of the free energy Fjq per unit 
length of the oo x N lattice for large N and for free and 
periodic boundary conditions. 

Let us first consider the case of an infinitely long strip 
of width N with free boundary conditions. For N even, 
the results of |^| show that 
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Fi (7V+l) + log(l + v / 2). 
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The expansion Q then gives 
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-£(JV+l) + log(l + >/2)-|™ + ... (5) 



Similar calculations for N odd lead to a formula in terms 
ofZ Q i(N + l), and yield 
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The analogous results for the periodic case, i.e. an infi- 
nite cylinder of perimeter N, read 
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The free energy per unit length of an infinitely long 
strip of finite width N at criticality has the finite-size 
scaling form |(J 



F = fbulkN + / sur f + — 



(9) 



where /bulk and / SU rf are, respectively, the bulk and the 
surface (boundary) free energy densities, and A is a con- 
stant. Though the free energy densities /bulk and / SU rf 
are not universal, the constant A is universal. The value 
of A is related to the conformal anomaly c of the un- 
derlying conformal theory, and depends on the boundary 
conditions in the transversal direction. These two depen- 
dencies combine into a function of the effective central 
charge c cff = c - 24/i min , 



A = -— c, 



eff 



, M 1 ' ' K h min -—j on a strip, (10) 

7T / C \ 

A = -— c ff = An [h m i n - —J on a cylinder. (11) 

The number h m i n is the (chiral) conformal weight of the 
operator with the smallest scaling dimension present in 



the spectrum of the Hamiltonian with the given bound- 
ary conditions (for the cylinder, we assumed that this 
operator is scalar, h m - ln = h m - ln ). 

In a unitary theory, one has h min — on a cylinder 
(with periodic condition) and on a strip with identical 
left and right boundary conditions, and h mm > other- 
wise. In a non-unitary theory, like the conformal theory 
discussed here, there is no restriction on h m i n . 

The finite-size corrections computed above have all the 
form @, with the effective central charge depending on 
the parity of N, see ifTUjl and <(TT|) . We will show that 
indeed the effective central charge, and not the central 
charge itself, depends on the parity of N, because value 
of ftmin does, due to the fact that changing the parity of 
N in effect changes the boundary condition. 

To understand this peculiarity of the dimer model, we 
consider, first on the strip then on a cylinder, the map- 
ping of the dimer model to the spanning tree model |l6j | 
and, equivalently, the Abelian sandpile model 

Dimers on a strip. - Let us consider first the dimcr 
model on the rectangular lattice C of size M x N with 
free boundary conditions. Since we are interested in the 
limit M — > oo, the parity of M will not matter here. For 
simplicity, we take M odd, and discuss successively the 
cases N odd and ./V even. 

For M and N both odd, the bijection between dimer 
coverings on C with one corner removed and spanning 
trees on the odd sublattice Q C C is well-known [3, LHj . 

A dimer containing a site of Q, in blue in Figure 1, 
can be represented as an arrow directed along the dimer 
from this site to the nearest neighbour site of Q. It is 
easy to prove that the resulting set of arrows generates 
a uniquely defined spanning tree, rooted at the corner 
which had been removed from C (see Figure 1). Since 
the dimers which do not contain a site of Q are com- 
pletely fixed by the others, one has a one-to-one corre- 
spondence between dimer coverings on C minus a corner 
and spanning trees on Q . This allows one to express the 
number of dimer configurations by the Kirchhoff theo- 
rem as Z — det A where A is the Laplacian matrix for 
spanning tree on Q. As shown in [Tt| . spanning trees on 
Q, rooted at a corner, are in bijection with the configu- 
rations of the Abelian sandpile model (ASM) on Q, with 
closed boundary conditions on the four boundaries, the 
only sink (dissipative) site being the root of the trees. 

When M — > oo, the lattice Q becomes an infinitely 
long strip of width N and, in the ASM language, closed 
boundary conditions on the two vertical sides. Many 
independent arguments and explicit calculations all con- 
verge to a central charge c = — 2 for the ASM on a square 
lattice 0, 0, . On the other hand, the spectrum of 
the ASM Hamiltonian on a slice of the strip with closed 
boundary conditions at the two ends has been computed 
in 18] . It is given by a single representation 1Z, reducible 
but indecomposable, of the chiral algebra of a rational 
logarithmic conformal field theory with c = — 2 |20|. The 
representation has two ground-states, the identity opera- 
tor and its logarithmic partner, both of conformal weight 
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FIG. 1: (Color) Mapping of a dimer covering to a spanning 
tree on the odd sublattice, for anMx7V = 5x7 lattice (top), 
and a 5 x 8 lattice (bottom). In both cases, the solid dots 
represent the sites of the odd sublattice Q, and the open dots 
are the roots of the trees. 



0, so that /i m in = 0. The effective central charge in this 
sector is therefore c c g = — 2, and the general formula (|1C)|) 
reproduces the finite-size corrections (0. 

When M is odd and N is even, dimer coverings exist 
on L without the need to remove a corner. In this case, 
the above construction leads to a set of spanning trees 
on the odd sublattice Q, where certain arrows may point 
out of the lattice from the right vertical side (Figure 1). 
Viewing this vertical boundary of Q as roots for the span- 
ning trees, we see that dimer coverings on C map onto 
spanning trees on Q which can grow from any site of the 
vertical side. In turn, the spanning trees map onto the 
ASM configurations with one vertical open, dissipative 
boundary, and the three other closed. 

In the limit M — > oo, the lattice becomes an infinite 
strip with open and closed boundary conditions on the 
two sides. In this case, the results of [3 show that 
the ground-state of the Hamiltonian with such bound- 
ary conditions is a primary field of conformal weight 
hmia = —1/8. With c = — 2, this yields c c fr = 1 and 
again the formula l|10|) correctly gives the result (J5J . 

Let us note that the bijection between the dimer cov- 
erings and the spanning trees holds if we use the even 
sublattice instead of the odd one. The boundary con- 
ditions however change. If N is odd, the vertical sides 
(and the horizontal ones as well for M odd) become open 
rather than closed. The spectrum of the corresponding 
Hamiltonian change, with a non-degenerate ground-state 
being the identity operator ^8| . Thus the value h m { n = 
remains. If N is even, the left and right boundaries, 
previously closed and open respectively, become open 
and closed, so that the Hamiltonian remains the same, 
hmia = —1/8. 

The equivalence of the odd and even sublattices is a 
duality property. The two sublattices are dual to each 
other, and the spanning trees on Geven are dual to those 
on Godd- It is not difficult to check that open and closed 
boundary conditions are exchanged under duality. 

Thus the leading finite-size corrections for an infinitely 



long strip of width N agree with the prediction of a 
c = — 2 conformal field theory, provided one realizes that 
changing the parity of N genuinely changes the bound- 
ary conditions, an effect due to the strong non locality 
of the dimer model. The change of boundary conditions 
is not apparent in the dimer model itself, but is manifest 
when one maps it onto the spanning tree model or the 
sandpile model. 

The primary field with h = — 1/8 can be further tested 
|l8j. If, on a rectangular lattice M x N, we remove n 
sites from the lower boundary, the height will be M or 
M—l and will therefore take on the two parities. This has 
the effect of changing the boundary condition along the 
lower boundary, from closed to open, and back to closed 
(or vice- versa). We have checked that the universal part 
of the ratio of partition functions, after and before the re- 
moval, and expected to be equal to (<j>h(0)<f>h(n)) ~ n~ 2h 
where <ph is the boundary field that changes a bound- 
ary condition from open to closed, is indeed asymptoti- 
cally equal to n 1 / 4 , in the limit M,N — > oo, supporting 
ft = -1/8. 

Dimers on a cylinder. - We consider here an M x N 
rectangular lattice C with periodic boundary condition 
in the horizontal direction, so that £ is a cylinder of 
perimeter N and height M . As before, we will eventu- 
ally take M to infinity, which makes its parity irrelevant. 
We choose M even. According to the discussion of the 
previous section, the top boundary of the cylinder then is 
subjected to open boundary conditions (in ASM terms) 
while the bottom one is closed. We separate the cases N 
odd and N even. If N is odd, we select the sublattice Q 
consisting of those sites of C having odd-odd coordinates. 
It is easy to see that two columns of Q will contain sites 
which are neighbors in Q and in C (connected by hori- 
zontal bonds). Therefore a dimer may contain zero, one 
or two sites of Q. The dimers containing no site of Q 
are completely fixed by the others and play no role. For 
the others, we do the same construction as before. If a 
dimer touches one site of G, we draw an arrow directed 
along the dimer from that site to the nearest neighboring 
site of Q . However, for a dimer containing two sites of G, 
the two arrows would point from either site to the other, 
ruining the spanning tree picture. It can nevertheless be 
restored in the following way. 

Instead of seeing the two arrows as pointing from one 
site to its neighbor, we say that they point towards roots 
inserted between the neighbor sites, thus replacing the ar- 
rows •— ► «— • by •— »o CM— •• This in effect amounts to 



cut the cylinder along the bonds of C which connect sites 
of Q, unwrapping it into a strip, and to add columns of 
roots on the left and on the right side of the strip. The 
new arrow configurations define spanning trees, rooted 
anywhere on the left and right boundaries. So dimer cov- 
erings on the original cylinder are mapped to spanning 
trees on a strip, with open top and closed bottom hori- 
zontal boundaries (we chose M even), and open vertical 
boundaries. 

When M goes to infinity, the lattice becomes an infi- 
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nite strip with open boundary condition on either side. 
As mentioned above, the ground-state of the Hamilto- 
nian is the identity, of weight h m in = 0, leading to an 
effective central charge c c fr = —2. The general formula 
(fTUfl for the strip gives the correct result iJSJ. 

This is a very unusual situation. Although the dimer 
model is originally defined on a cylinder, it shows the 
finite-size corrections expected on a strip, and must really 
be viewed as a model on a strip. 

For N even, the problem of having two arrows pointing 
from and to neighbor sites does not arise, however the ar- 
rows one obtains do not define spanning trees, but rather 
a combination of loops wrapped around the cylinder and 
tree branches attached to the loops. Each loop has two 
possible orientations. The one-to-one correspondence be- 
tween the oriented loops combined with tree branches 
from one side and dimer configurations from the other 
side can be established as above. The enumeration of 
the loop-tree configurations needs a generalization of the 
Kirchhoff theorem, Z — det A where A is the discrete 
Laplacian with antiperiodic boundary conditions. In the 
continuous limit, this leads to the free theory of antiperi- 
odic Grassmann fields which, in turn, gives h m i n = — 1/8 
|2lj . From the general formula for the cylinder, we 
see that the finite-size correction (Q) is again consistent 
with c = —2. 

The above interpretation of the peculiarities of the 
dimer model shed a new light on old calculations by Fer- 
dinand of the part ition function of the dimer model on a 
M x N torus [23. When M and N are both even, the 



universal part of the partition function equals 

_ ej + el + g| 

in the limit M, N — > oo with fixed ratio q = exp (— 2tt^-). 
This is exactly equal to the partition function Z c= -2{q) = 
Ix-i/sl 2 + 2|xo + Xi\ 2 + IX3/ 8 |l°f the c = -2 rational 
conformal theory developed in [2(J , confirming the value 
of c = —2 (and the value h m i n = — 1/8 found for the even 
cylinder). In case one of M or N is odd, the partition 
functions are 

Zodd,even = 7T~ (a/?)) Z even o( id — — (13) 

zr] 2ry 

These are cylinder partition functions of the same c = — 2 
conformal theory |l8j . in agreement with the view that 
a periodic dimension of odd size is actually not periodic, 
and correspondingly tori with one odd dimension are in 
fact cylinders. 
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